Intensity fluctuations in steady state superradiance 
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Alkaline-earth like atoms with ultra-narrow optical transitions enable superradiance in steady 
state. The emitted light promises to have an unprecedented stability with a linewidth as narrow as 
a few millihertz. In order to evaluate the potential usefulness of this light source as an ultrastable 
oscillator in clock and precision metrology applications it is crucial to understand the noise properties 
of this device. In this paper we present a detailed analysis of the intensity fluctuations by means 
of Monte-Carlo simulations and semi-classical approximations. We find that the light exhibits 
bunching below threshold, is to a good approximation coherent in the superradiant regime, and is 
chaotic above the second threshold. 

PACS numbers: 42.50.Fx, 37.30.+i, 42.50.Pq, 42.50.Ct, 42.55.Ah, 42.50.Lc 



I. INTRODUCTION 

Dicke superradiance is a paradigmatic collective effect 
in cavity quantum electrodynamics (QED). At the fun- 
damental, level superradiance is a quantum interference 
effect in which the probability amplitudes for the emis- 
sion of a photon from several atoms conspire to yield 
a collective light emission rate that is larger than for 
uncorrelated atoms. Due to its great conceptual sim- 
plicity and generality superradiance has been extensively 
studied both experimentally and theoretically. The noise 
properties of the light emitted in superradiance have been 
of particular interest. This is because the early stages 
of superradiance are often initiated by quantum fluctua- 
tions which are subsequently amplified by the collective 
emission process [l|, |2| . Superradiance can thus serve as a 
physical phenomenon that allows us to study the micro- 
scopic quantum fluctuations through their macroscopic 
consequences. Examples of this macroscopic manifesta- 
tion of the quantum fluctuations are the first passage 
time statistics of the superradiant pulse ^Q as well as 
the second order correlations of the field J61| . 

Nearly all realizations of Dicke superradiance have 
been in the pulsed regime. In contrast to such exper- 
iments we have recently proposed a system based on 
earth-alkaline-like atoms in which superradiance can be 
achieved in steady state [3, Q • The interest in that light 
source derives from the extremely narrow linewidth of the 
generated light. For experimentally realizable parame- 
ters linewidths in the millihertz range could potentially 
be realized. The light generated this way could thus serve 
as an ultra-stable local oscillator with a stability that is 
about two orders of magnitude better than the current 
state-of-the-art. At the core of this device are atoms 
with an ultra-narrow optical transition coupled to a high 
finesse cavity. The atoms collectively emit photons into 
the cavity mode and they are concurrently repumped to 
the excited state, providing a steady supply of energy. 
The collective decay of the atoms via the cavity mode es- 
tablishes a collective atomic dipole, which radiates much 
more strongly than independent atoms would. Depend- 
ing on the repumping rate, the system can also exhibit 



subradiance, or thermal light as would be the case for 
an ensemble of random radiators. Qualitatively similar 
behavior has been predicted for the overdamped many- 
atom micromaser [9|. 

Just as for pulsed Dicke superradiance, the higher or- 
der correlations of the light are non-trivial and a solid 
knowledge of them is crucial for a full understanding 
of the collective light generation mechanism as well as 
for potential applications. Correlations of the intensity 
can be used to quantitatively study these fluctutations 
and these are the subject of this paper. Some aspects 
of the noise properties of continuously pumped, collec- 
tively emitting systems have also been discussed recently 
in the context of collectively radiating low dimensional 
solid state systems [Ifl, \lM- Fo'" instance Temnov and 
Woggon studied the photon statistics deep in the subra- 
diant regime in [12J. 

The goal of this paper is to fully characterize the in- 
tensity correlations of the light generated by means of 
steady state superradiance. To this end we introduce a 
simplified model that captures all the essential aspects 
of the problem in section |lT] and we recall the basics of 
steady state superradiance. Section IIIII presents our re- 
sults on the Hanbury Brown-Twiss correlations of the 
generated light obtained for small atom numbers using 
quantum Monte-Carlo simulations and semiclassical ap- 
proximations in the limit of large atom numbers. 



II. MODEL 

Calculating higher order correlations of quantum fields 
is typically a hard problem. Analytic closed form solu- 
tions are known in only a few special cases and thus we 
have to rely mostly on numerical simulations. Because 
of the exponential scaling of the dimension of the Hilbert 
space of the system with the number of atoms we must 
restrict ourselves to the simplest model that still captures 
the essential physics that we are interested in. 

The core ingredients of a system exhibiting steady 
state superradiance are illustrated in Fig. [T] An ensem- 
ble of N two level atoms with excited state |e), ground 
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FIG. 1: (Color online) Schematic illustration of N two level 
atoms confined in a single mode cavity field indicated by the 
black dots. The atoms are being incoherently repumped. The 
output of the cavity field is monitored by two detectors Di 
and D2 with a variable time delay r between them. 



state I (7), and transition frequency Wq are collectively cou- 
pled to a high finesse cavity with resonance frequency uJc- 
The atoms are independently repumped from the ground 
state to the excited state in order for the atoms to be 
able to radiate photons continuously. 

The non-collective nature of the repumping is essential 
for two reasons. First, it is much easier to achieve experi- 
mentally than collective repumping, for instance by opti- 
cal pumping through an auxiliary excited state. Second, 
and somewhat paradoxically, it can balance the effects 
of other incoherent processes such as spontaneous emis- 
sion and dephasing that would normally drive the atoms 
into less collective states that do not exhibit superradi- 
ance. A p urely collective repumping as considered e.g. 
in |13l . [l4| cannot change the length of the atomic Bloch 
vector. Therefore, if the length decays due to dissipative 
processes, it cannot be restored and the superradiance 
stops. In contrast, the length of the atomic Bloch vector 
can grow in the case of non-collective pumping so that a 
collective atomic dipole can develop from an ensemble of 
completely independent atoms. 



A. Atom-field master equation 

Mathematically, the coupled atom cavity system can 
be described by the Hamiltonian 

H = hwaJz + hwca*a H {a* J^ + J+a) . (1) 

Here we have introduced an angular momentum rep- 
resentation for the atoms in the usual way by identifying 
the excited state with the spin up state of a fictitious 
spin 1/2 system and the ground state with the spin down 
state. The operators 
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are the z-component and ladder operators of the total an- 
gular momentum. In these equations cfz = l^) (^l ~ I5) (5I 
is a Pauli matrix pertaining to atom number j and 
(T_ = cr\. = \g){e\ is the spin-fiip operator for atom j. 
The operators a and a^ are the annihilation and creation 
operators for a photon in the cavity mode. For simplicity 
we assume that the coupling constant g of the atoms to 
the cavity is identical for all atoms. This could, in princi- 
ple, be achieved by trapping the atoms at the antinodes 
of the cavity mode. Less ideal spatial configurations of 
the atoms merely lead to a reduction of the effective num- 
ber of atoms by a factor of order one that clearly has no 
impact on the basic conclusions of this paper. 

Besides the coherent interaction of the atoms with the 
cavity field there are also dissipative processes due to the 
coupling of the atoms to field modes outside of the cavity 
and due to decay of the cavity fields. The decay of the 
cavity can be accounted for with the usual Born-Markov 
master equation for the reduced density matrix for atoms 
and cavity field, 

^ = ^ [p, 7^] + /:eav [p] + ^pump [p] , (4) 

where the Liouvillian for the cavity decay with intensity 
decay rate k is 

'Ccavi/o] = --{a)ap + pa) a - 2apa)) . (5) 

The repumping of the atoms with pump rate w is de- 
scribed by the Liouvillian 
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pump 
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(6) 

We are assuming that the spontaneous emission of the 
atoms into free space with decay rate 7 can be ne- 
glected. In general this assumption requires C ^ 1 where 
C = g"^ /{n'y) is the single atom cooperativity parameter, 
see Eq. ^ below. Note however that in the superradiant 
regime where the decay through the cavity is collectively 
enhanced it is found that the much less stringent condi- 
tion NC 3> 1 is sufficient for this approximation to be 
justified. 

In order to calculate correlation functions of the gen- 
erated light field we simulate the dynamics of the system 
subject to the master equati on ffl using the Monte-Carlo 
wavefunction technique [15l4l7l |. In that technique the 
evolution of the system is represented by an ensemble of 
stochastic wavefunction trajectories {|'0(^))} where each 
trajectory \i){t)) is a representative evolution of the sys- 
tem. 



B. Adiabatic elimination of field in "bad cavity" 
limit 

The Hamiltonian Eq. ([T]) is suitable for the Monte- 
Carlo simulations because it directly grants access to the 
field correlations we are interested in. For analytical cal- 
culations it is desirable to further simplify the problem by 
exploiting the fact that the cavity field decays so much 
faster than the atomic coherence. Adiabatically elimi- 
nating the light field yields the effective superradiance 
master equation |6|], 



tion values of the field amplitude, 



dp_ _ _y_c 

dt~ 2 



J+J-P + pJ+J- - 2J-pJ^ 
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Here, the collective decay rate of the atoms is given by 

r, = C7 = gV'* ■ (8) 

The condition for the validity of the adiabatic elimina- 
tion of the cavity field is that the cavity field relaxes 
much faster as the atoms. Using that the fastest atomic 
relaxation rates are obtained in the superradiant regime 
and that they are of order NC^ we find the quantitative 
condition 



K » NC-1 . 



(9) 



In this "bad cavity" limit the role of the cavity mode 
reduces to providing a collective decay channel for the 
atoms. The simplification brought about by the elimina- 
tion of the field is two-fold. It allows us to deal with the 
atomic degrees of freedom only and all parameters of the 
full coupled atom-cavity system have collapsed into just 
one characteristic parameter, w/Tc- 



C. Hanbury-Brown-Twiss signal 

Fluctuations of the light intensity can be characterized 
experimentally by a Hanbury-Brown-Twiss-like setup as 
illustrated in Fig. [1] [ll, [ll]- The light passes through a 
50/50 beam splitter and the intensities in each output- 
port of the beam-splitter is detected with a photo-diode. 
The photo-currents of each detector are then multiplied 
and integrated. A variable delay r can be imparted on 
one of the outputs of the beam splitter in order to mea- 
sure correlations of the field at different times. Using 
such a setup it is possible to measure the joint probability 
P2{t,t + T)AtAT to detect a photon both in a time inter- 
val At at t and in a time interval At at i -I- r. According 
to the theory of photo-detection [T^ [20| this probability 
can be calculated in terms of normally ordered expecta- 
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Pl{t)AtPl{t + T)AT 

{a) {t)a^ {t + T)a{t - 
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(at(i)a(t))(at(i-hT)a(i-hr)) 



(10) 



We have normalized the joint probability to the sin- 
gle time probabilities for photon detection, Pi{t)At and 
Pi{t + At) At. In steady state, which is the case of 
primary interest here, g^'^'{t,T) does not depend on t. 
For notational convenience we drop the variable t from 
(/'^'(t, t), i.e. we simply write 5'^' (r). In writing Eg. (|10p 
we have also made use of the result from the input- 
output theory for cavities that normally ordered corre- 
lation functions outside of a cavity are equal to the nor- 
mally ordered correlation functions of the intra-cavity 
field, provided that the input ports of the cavity are in 
vacuum states [2l|, ^^ . 

The second order correlation function at zero time de- 
lay, 5^^^(0), is related to the fiuctuations AP = {P) — 
{I)^ of the out-coupled photon flux /. These fluctuations 
AP can be used to characterize the instability of the in- 
tensity of the out-coupled beam because / is proportional 
to the beam intensity. Typical photo detectors have a de- 
tection bandwidth B that is extremely large compared to 
the cavity bandwidth and in that case the fluctuations of 
the intensity are given by 



AP 



/2(.g(2)(0)_l) + B/ 



(11) 



where I — (/). For coherent light the arrival of pho- 
tons at the detectors are a Poisson process in which 
the arrival times are completely random with mean rate 
of arrival /. For such light AP = BI and we have 
5^^-'(0) — 1. Light with larger intensity fluctuations is 
called super-Poissonian and light with smaller intensity 
fiuctuations is called sub-Poissonian. Super-Poissonian 
light has g^'^> (0) > 1 and photons arrive in bunches, while 
sub-Poissonian light has 5*^^^ < 1 and photons are anti- 
bunched, i.e. they arrive more regularly than pedicted 
by Poisson statistics. 

Our numerical wave- function Monte-Carlo simulations 
grant us access to correlation functions in two ways. 
First, we can calculate expectation values of system ob- 
servables O in the usual way by calculating (^jOIV') and 
averaging over the ensemble of trajectories. Alternatively 
we can extract the correlation functions by an analysis of 
the decay times of the system that very closely mirrors 
an actual experimental procedure. This latter method is 
easier to implement for non-zero time delays. All results 
on the intensity correlations presented here were calcu- 
lated this way due to the greater fiexibility of having 
access to zero and non-zero delays. 

The procedure by which we calculate g^ \t) from the 
cavity decay times is illustrated in Fig. [2] Part a) of that 
figure shows the evolution of the mean photon number 
inside the cavity for an example trajectory. Cavity decay 
events are indicated by the black ticks at times ti. For 
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FIG. 2: (Color online) Time evolution of the intracavity in- 
tensity during a single Monte-Carlo trajectory (a) for A'^ = 10 
atoms and w — bTc- The black ticks indicate cavity decay 
events. Panel (b) shows the binning of the decay times sub- 
sequent to ti leading to the histogram nij(At) used in the 
calculation of (;'^'(r) (see text for explanation). 



each photon emission event we calculate the Histogram 



'*j 



(Ai) = #photons in {U + jAt,U,+{j + l)Ai]. By 
averaging these histograms over all i we find nj(Ai) — 

'^phot Si "■i.i(^0: where riphot is the total number of pho- 
tons emitted. This histogram is closely related to the con- 
ditional probability P{t+jAt\t) to find a second photon 
at time t + jAt provided that a first photon has been 
detected at time t, i.e. nj{At) = AtP{t + jAt\t). Us- 
ing the relation between that conditional probability and 
the joint probability, P2{t,t + jAt) = P{t + jAt\t)Pi{t), 
we can then find g^'^'{jAt) on a grid with spacing At 
according to 



9^'Hj^t) 



^phot/#bins ' 



(12) 



where nphot/^T^bins = Pi{t)At is the mean number of 
photons per bin. The choice of the bin width At is a 
trade off between resolution and statistical fiuctuations 
and it has to be chosen differently for different simula- 
tion parameters. It must be small enough to resolve the 
dynamics of the system. Once that constraint is satis- 
fied it should be as long as possible in order to yield the 
smallest fluctuations. 



III. RESULTS 

Generally, as explained in detail in [8|, three regimes 
of light emission can be distinguished depending on the 
repump strength. If the repunip rate is smaller than the 
effective atomic decay rate, w <C Fc, the atoms evolve 
into a dark state in which the emission of photons is 
strongly suppressed despite nearly half the atoms being 
in the excited state. In the intermediate regime, Tc < 
w < NTc, the atoms emit light in a superradiant fashion 




FIG. 3: (Color online) Second order intensity correlation 
g^ '{0) as a function of the repump rate. The green symbols 
show the Monte-Carlo results including the statistical errors 
for TV = 10 atoms. The blue solid line shows the analytical 
result Eq. (|17|) for A^ = 10 atoms, the purple dashed line 
is for A^ — 100 atoms, and the yellow dash-dotted line is for 
iV = 1000 atoms. The gray dashed lines at g'^'(O) = 2 and 
g (0) — 1 are for orientation. 



and in the strong pumping limit, w 3> NTc, nearly all 
atoms are in the excited state and they emit chaotic light 
like an ensemble of thermally excited atoms. 



A. Monte-Carlo results for g*^'(0) 

First we consider g'-^'{0) shown in Fig. [3] for A^ = 10 
atoms. The error bars in that figure are estimates of 
the statistical uncertainty obtained by treating the his- 
tograms nij{At) for different i as independent of each 
other. 

In the weak pumping limit the light exhibits strongly 
super-Poissonian fluctuations indicating photon bunch- 
ing. This bunching effect can easily be understood in the 
extreme limit w/Tc -^ ^, U^ . In that limit the atoms 
are optically pumped into collective dark states [23, [23] 
I J = 1,M = -1) and |J = 0,M = 0) '2^. From 
\J = 0,M = 0) the atoms can only be pumped to | J = 
1, M = 1) from which they relax to \J — l,M ^ —1) by 
rapidly emitting a pair of photons in a cascade within a 
time of order T~^. 

In the superradiant regime for intermediate pumping 
.g'-^'(O) is reduced to a value about half way between 1 
and 2. It reaches a minimum at the superradiant emission 
maximum, w — NTc/2. The minimum value depends on 
N with larger N yielding smaller g'^^^(O). 

In the limit of very strong repumping the atoms behave 
like a thermal ensemble of atoms. Therefore we have 
g'-^-' (0) = 2(1 — 1/A^) in that limit in very good agreement 
with the Monte-Carlo simulation results. 



B. Semi-classical results for g^^'{0) 

In applications, the number of atoms will typically be 
much larger than 1. That makes it easier to achieve the 
collective strong coupling regime where NC ^ 1. Fur- 
thermore, the collective nature of the light emission is 
more apparent in that limit. 

Unfortunately, the Monte-Carlo simulations on which 
the above results were based, cannot be easily imple- 
mented for the study of large atom numbers. The rea- 
son for this is that the size of the system Hilbert space, 
d = 2^, scales exponentially with the number of atoms. 
In this section, in order to bypass these difficulties, we 
exploit the possibility for a semiclassical approximation 
that precisely derives from the largeness of N. The key 
idea is that in a macroscopic ensemble of atoms the cor- 
relations between n atoms can be expressed in terms of 
the correlations between n — 1 atoms plus an error term. 
Ordinarily, the error terms become smaller as the "clus- 
ter size" n is increased. The approximate treatment that 
we employ here assumes that correlations between more 
than two atoms can be completely expressed in terms of 
pair correlations and single atom quantities. 

The semi-classical calculation involves two non-trivial 
steps. First we have to find the correlations between 
the spins of different atoms. Second we have to find an 
expression for g^"^' (0) in terms of the atomic correlations. 
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FIG. 4: (Color online) (a) Inversion (ctz ), (b) spin-spin cor- 
and (c) spin-spin correlation {ai a\ ) - 



relation {a^ a_ ), 

(ai )^ as a function of pump strength for A'' = 10 (blue solid 



line), A'^ = 100 (purple dashed line), and A^ 
(yellow dash-dotted line). 
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1. Steady state solutions for pair correlations 

The symmetry of the expectation values with respect 
to particle exchange greatly reduces the number of ex- 
pectation values that have to be considered. We have for 
instance {a^^'^a^l'') == ((7^^'ai^^) for ah i ^ j. Up to the 
level of pair correlations, all observables that we are in- 
terested in can be expressed in terms of (o-^ ) , {a\_ (t_ ) , 
and (ai^Vi'^). 

The equations of motion for these expectation values 
can be found from the master equation ([7|), 



dt 



-(^ + r,)((aW)-do) (13) 



2Fe(iV-l)(4^Vi'') 



where da — {w ~ Vc)/('w + Fc), 

J/-(1)'(2)\ 



dt 



-(z«-fF,)(<7WaL')) 

+ F.(^-2)(aWal^)ai^)), 



(14) 



dt 



-2(^. + ^,)((<7Wai2))-do(aW)) (15) 
+ 4F.((a«ai^)) 

-(iV-2)(^«4^)a(^)) 



We have checked that the third order expectation val- 
ues {(Tz (t!^ al ) can be factorized according to 

(a«ai^)^i^))«(aW)(a(.^)ai^)) 

to a very good approximation by evaluating them in our 
Monte-Carlo simulations and by approximating them in 
terms of lower order cumulants in the adiabatic approx- 
imation. By factorizing this way we obtain a closed set 
of equations. The steady state expectation values are 
obtained by setting the time derivatives equal to zero. 
The resulting algebraic equations can be solved analyti- 
cally for {(Tz ), ((T4_ (T_ ), and {(Jz oz ), leading to rel- 
atively complicated expressions that we reproduce in the 
appendix for completeness. Plots of the steady state ex- 
pectation values are given in Fig. U] for different atom 
numbers. For large atom numbers the inversion (ctz) is 
essentially zero below threshold for collective emission, 
increases linearly with w in the superradiant regime, and 
saturates with all atoms in the excited state in the strong 



pumping regime. In the superradiant regime, the spin- 
spin correlation {a\. a_ ) is approximately described by 
an inverted parabola with zeros at the thresholds w — Tc 
and w = NTc and a peak value of 1/8 at w = NTc/'2,. 
In the strong pumping regime these correlations are de- 
stroyed by the repumping. The spin-spin correlations 

{c^z cFz ) approximately factorize as (iri '(Jz ) ~ {(^z )^ 
in the limit of large atom numbers. Note that the equa- 
tions for {ai ' ) and (a)^ a_) close if that factorization 
is made and consequently much simpler approximate ex- 
pressions can be obtained for w > Fc in the limit iV — ^ oo 
as pointed out in Q. However, since we would like to 
compare the semi-classical theory with the Monte-Carlo 
results that were obtained for only relatively small atom 
numbers, we have to use the more complicated expres- 
sions discussed here. 



2. Expression of g^ \Q) in terms of atomic operators in the 
had cavity limit 



In the bad cavity limit discussed in subsection lllBl the 
cavity field is slaved to the atomic dipoles. 



in 



(16) 



This means that we can calculate correlation functions 
of the field if we know atomic correlation functions, for 
instance 

J2).^^ _ {J+J+J-J-) 



g^'Ko) 



{J+J-) 



^{i)c^U)^ik)c^il)\ 






(17) 



The atomic expectation values can be expressed in terms 
of the above expectation values. 



N 



J2 (afd^) = N{{a(P) + l)/2 + N{N l)(aV^'ai^') , 

i J — 1 

(18) 
and 

{J+J+J-J-) = N{N-l) (19) 

x(2(iV-2)((aW) + l)(aV^)aL')) 

+ (l + (aWa(^))+2(aW)2)/2 
+ (iV-2)(Ar-3)(a^''ai'')2) . 

In order to arrive at this last result we have factorized 
expectation values for four different atoms according to 
(<5'4_ (T4_ (T_ CT_ ) w ((T^^(7_ )^. Expectation values in 
which at least two indices are identical involve at most 
three different atoms and, with the factorization dis- 
cussed earlier, they reduce to the known pair correlations 
and single atom expectation values. 



3. Comparison with Monte-Carlo results 

The semi-classical results for ^'■^■'(O) are also shown in 
Fig. [3] The semi-classical curve agrees very well with 
the Monte-Carlo results for w > Tc- Below that thresh- 
old the semi-classical expression yields unphysical values. 
The disagreement below threshold is not surprising be- 
cause the atoms are in a very highly correlated state in 
that regime and these correlations cannot be captured by 
taking into account only pair-wise correlations. 

The good agreement between semi-classical and 
Monte-Carlo results for w > Tc allows us to use the 
semi-classical expression to extrapolate to very large 
atom numbers. We find that for large atom numbers 
the field exhibits nearly coherent counting statistics, i.e. 
g^^\0) w 1, in the superradiant regime, Tc <^ w < NTc, 
and it has the counting statistics of chaotic light in the 
strong pumping regime. 



C. Monte-Carlo results for g (t) 

The correlations between photons arriving at a pho- 
todetector characterized by g^'^\0) only persist for a cer- 
tain amount of time. To study the decay of these inten- 
sity correlations we show the second order correlation for 
non-zero delay, r 7^ 0, for 10 atoms in Fig. [Sj 

The strong bunching peak in the subradiant regime 
decays on a time scale of order Fj^ ^^. After that an 
anti-correlation dip develops for a period of order ~ w~^ 
because the repumping has to take the system out of the 
dark states before another photon can be emitted. 

In the superradiant regime the bunching is much 
weaker. The bunching peak also disappears on a much 
shorter time scale of order ~ l/(iVFc). The subsequent 
anti-correlation is a much weaker effect that disappears 
as 1/A^ in the limit of large atom numbers. For times 
T ^ w^^, (A^Fc)"^ the light intensities are uncorrelated. 
This is remarkable because the amplitude of the field is 
coherent for much longer times of order F^^. Intensity 
correlations decay on a time scale set by the collective de- 
cay of the system while the decay of the amplitude corre- 
lations occurs on the time scale set by the single particle 
decay. From the perspective of potential applications of 
this light source as an ultra-stable local oscillator this 
means that the field can be considered as coherent to a 
very good approximation. 

As pointed out above the atoms behave like a thermal 
ensemble in the strong pumping regime. The two-time 
correlation function agrees well with the result for ther- 
mal light emitted by a large number of atoms [25| , 



g^'>{t)^i + \g'''{r) 



(1) 



l + e 2^/" 



(20) 



The slight discrepancy between this formula and the nu- 
merical results is likely due to the relatively small number 
of atoms considered while Eg. ([20]) is derived in the limit 

A^^ 00. 
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Such excess fluctuations could adversely affect the utility 
of light sources based on steady state superradiance as a 
stable frequency reference, one of the main motivations 
for studying this system in the first place. For such ap- 
plications it is crucial that the long coherence time and 
collectively enhanced brightness demonstrated previously 
[7| are paired with small intensity fluctuations. 

In the subradiant regime the emitted light exhibits 
strong bunching and super-Poissonian intensity fluctu- 
ations. As has been pointed out previously by Temnov 
and Woggon [T^ this effect could be useful in identifying 
and analyzing the subradiant regime in experiment. The 
well understood thermal character of the atomic ensem- 
ble in the strong pumping regime serves as a valuable 
benchmark for the validity of our theoretical treatment. 

In future research we plan to systematically investigate 
the cross-over from steady state superradiance to a laser. 
In the extreme limits of this cross-over the system is dom- 
inated by purely atomic collective enhancement on the 
one hand and by purely photonic collective enhancement 
through stimulated emission on the other hand. The in- 
termediate regime between the two were both collective 
enhancement due to stimulated emission and an atomic 
collective state are equally important is very intriguing 
from a fundamental point of view. 
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FIG. 5: (Color online) Second order correlation g (t) for 
N = 10 atoms (a) for the subradiant regime for weak pumping 
w — 0.25rc, (b) for the superradiant regime with w = S.OFc, 
and (c) for the strong pumping regime with w = lOOFc. Note 
the different scales on the time axis for each figure. The purple 
line in (c) is the thermal light result. 



Acknowledgments 

We are grateful for stimulating discussions with J. K. 
Thompson, Jun Ye, and J. Cooper. This work has been 
supported by NSF and ARO. 



Appendix A: Steady-state atomic pair correlations 

In this appendix we summarize the analytical solution 
of the semi-classical equations Eqs. (|13m5p . 
In steady state, the atomic inversion is 



.«)=cio-^^^^^^(^iV^V 



(Al) 



Inserting that into Eq. (|15|) we end up with a linear equa- 
tion for {&z (Jz )■ Solving that equation for {oz Oz) 
yields 



IV. CONCLUSION 

The key result of this paper is that the light emitted 
in steady state superradiance is second order coherent 
in the limit of large atom numbers. This result is sig- 
nificant because it establishes that the coherence prop- 
erties of the emitted light are closer to those of a laser 
than those of light generated in ordinary pulsed superra- 
diance. In contrast, light generated in pulsed superradi- 
ance would have super-Poissonian intensity fluctuations. 



2Tr 



((u; + Fc)(l-do(2iV-3)) 
+ 2(7V-l)(iV-2)r,((T^^^ai2))) (A2) 



(^. ^. )-d0 ' (^ + r,)2 



Inserting Eq. (IA1[) and Eq. (IA2[) into the remaining 

equation for (ct^ (T„ ) leads to a quadratic equation. One 
of the solutions must be discarded because it violates 
|((7!|_ (7_ )c| < 1 for certain repump rates and hence is 
unphysical. The physically acceptable solution is 



(^«^i^)) 



w + Tc 



A{N - 1){N - 2)wTl 



w^ + {2-{N- 2)do)wTc + (iV - 1)(1 + do)Tl 



- yidoil + do){N -1){N- 2)wTl + (u;2 + (2 - (iV - 2)do)wT, + {N - 1)(1 + do)Tiy . (A3) 

I 



This solution can then be inserted in Eq. (jAl[) and Eq. 
[U to find {aP)c and (ai^^ai^^)c. 
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